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Abstract. We introduce a sub-Riemannian analogue of the Bence-Merriman-Osher algorithm [41] 
and show that it leads to weak solutions of the horizontal mean curvature flow of graphs over sub- 
Riemannian Carnot groups. The proof follows the nonlinear semi-group theory approach originally 
introduced by L. C. Evans [27] in the Euclidean setting and is based on new results on the relation 
between sub-Riemannian heat flows of characteristic functions of subgraphs and the horizontal 
mean curvature of the corresponding graphs. 

1. Introduction 

The study of mean curvature flow in the sub-Riemannian setting is still at a very early stage, 
with several key properties, such as existence, uniqueness and regularity, still unknown. The notion 
itself of motion by mean curvature is understood only in special cases, as in the evolution of graphs 
over groups. 

In this paper we give two contributions to this topic. First we establish in the Carnot group set- 
ting a formula relating mean curvature of a graph and the heat flow of the characteristic function of 
the corresponding sub-graph (see Lemma[3J). This formula seems to be new also in the Riemannian 
setting (see Corollary [1]) . Then, using such formula we prove convergence of an analogue of the 
Bence-Merrriman-Osher [41J algorithm in the Carnot group setting (see Theorem [1]), yielding a eas- 
ily implementable time-step approximation of the sub-Riemannian mean curvature flow of graphs 
over Carnot groups. This convergence can in turn be used to define a notion of mean curvature 
flow also across characteristic points. 

1.1. Carnot groups. A Lie group G is called a homogeneous stratified Lie group if its Lie algebra 
admits a stratification 

g = V 1 © • • • V r with [V\ V j ] = V i+j and [V\ V r ] = 0. 

We will let n denote the topological dimension of G. Given a positive definite bilinear form go on 
V 1 we call the pair (G,go) a Carnot group and the corresponding left invariant metric go a sub- 
Riemannian metric. For each X 6 V' 1 let d(X) = i be the degree of X and let Q = Y2l=i * dim(Vi) 
the homogeneous dimension of G. Let us consider an orthonormal basis {^i, ....,X n } of Q fitted 
to the stratification, i.e. such that 

(1.1) the horizontal layer Ho ■= V 1 of G is generated by X±, . . . ,X m , 
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while {Xi}d(Xi)=k generates V k . These assumptions allow to define homogeneous dilations and an 
homogeneous norm on the tangent space at every point. If v E G is expressed in term of the element 
of the basis, 

n 
i=l 

In this setting we call horizontal gradient of a function u : G — > R the vector 

(1.2) V « = (Xiu, • • • X m u), 

and we say that u is in C\ if its horizontal gradient exists and it is continuous. The graph of a C\ 
function u : G — > R can be considered a surface in G x R. The product GxM can be endowed with 
the structure of Lie group and in particular with a Carnot group structure by setting X n +\ = d n +i, 
and d{n + 1) = 1, 

(1.3) the horizontal layer Hq = Hq x R of G x R is generated by (X\, ■ ■ ■ , X m ,X n+ i). 

The metric go is extended to a metric go in Hq by requesting that X n+ \ is orthogonal to Q and has 
unit norm. We call horizontal normal of the graph of u the projection on the horizontal plane Hq 
of the Euclidean normal and set 

C-Vo«,l) 



^i + |V(H 2 

Let us explicitly note that such graphs do not contain any characteristic points, i.e. points where the 
horizontal space is contained in the tangent space. In the literature several equivalent definitions 
of the horizontal mean curvature ho have been proposed at non characteristic points of a C 2 
surface M. To quote a few: ho can be defined in terms of the first variation of the area functional 
|191 [3"T| [T3l l4"6l l4"9l |4"3"1 111] , as horizontal divergence of the horizontal unit normal or as limit of 
the mean curvatures h e of suitable Riemannian approximating metrics g e |11| . If the surface is not 
regular, the notion of curvature can be expressed in the viscosity sense (we refer to [5], [6], |51| . 
|52j . [39] . [2], [ID], [9] for viscosity solutions of PDE in the sub-Riemannian setting). 
The curvature ho of the graph of a function u can be written as 

The horizontal mean curvature flow of a graph over a Carnot group G is the flow t — > Mt := 
{(x,u(x,t))\x G G} C G x R in which each point in the evolving manifold moves along the horizontal 
normal, with speed given by the horizontal mean curvature. The evolution of a family of graphs of 
functions t — > u(-,t) is then characterized by the following equation 

(1.5) d t u = -Au, where Au = - £ ySjj - - ^-^-^ jXjXjU, 

Here we will consider viscosity solution of this equation: 

Definition 1. A continuous function u : G x (0, oo) — > R is a weak sub-solution (resp. a super- 
solution) of equation (|1.5p if for every (x, t) S G x R and every smooth cp : G x (0, oo) — > R such 
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that u — <f> has a maximum (resp. a minimum 




d t (j)(x,t) < (resp. >) ^ (s,, 



) 



XiX^{x,t). 



Solutions are functions which are simultaneously super-solutions and a sub-solutions. 

Existence and uniqueness of viscosity solutions to this equation, attaining an assigned initial 
condition / has been established in [9] (see also [25] for a probabilistic interpretation of the flow) . 
One of the goals of this paper is to provide a discrete approximation of this motion, called diffusion 
driven motion by mean curvature. 

1.2. Diffusion driven motion by mean curvature. In the Euclidean setting, the motion by 
mean curvature can be obtained through an algorithm introduced by Merriman, Bence, Osher [41] . 
which relates the mean curvature flow and the heat flow. 
The algorithm is organized in two steps: 

• given A G K and a function /, the characteristic function xs x °f eac h of its A— level sets 
S\ = {x : A < f(x)} is diffused for a time t, via the subriemannian heat flow, giving rise to 
a smooth function w\(x, t). 

• at time t one defines a new function H(t)f by requiring that its A— level sets are {x : 



This procedure leads to the definition of a two step algorithm, which given a function / allows 
to define a new function H(t)f. The mean curvature flow t —> u(-,t) with initial data /, can be 
recovered applying iteratively this two step algorithm for shorter and shorter time intervals (or 
equivalently through iteratively applying the operator H). 

The convergence of the Euclidean version of this scheme has been proved by Evans [27] and 
Barles, Georgelin [3]. In particular Evans gave a proof based on nonlinear semi-group theory and 
on a pointwise study of heat flow of sets in terms of curvature of their boundaries. We also refer 
to Ishii [34], Ishii, Pires, Souganidis [35], Leoni [38] and Chambolle, Novaga [12] for extensions of 
this algorithm to more general setting. Let us mention that other convergent schemes for nonlinear 
parabolic equations have been proposed in [18], [20], [21], [22], [50] and [26] . 

1.3. Diffusion driven motion by horizontal mean curvature in Carnot groups. One of 

our main results is an extension of the algorithm in [H] to the degenerate parabolic setting of sub- 
Riemannian Lie groups. A motivation for such extension comes from problems of visual perception 
and modeling of the visual cortex. Geometric models for the visual cortex as a contact structure go 
back to work by Hoffmann [33], Petitot and Tondut [15] and Citti, Sarti and Petitot jH]. Later, in 
[15j . Citti and Sarti endowed this contact structure with a sub-Riemannian metric, and proposed 
a model of perceptual completion based on two cortical mechanisms applied in sequence. This 
two-steps process lead to a sub-Riemannian diffusion driven motion by horizontal mean curvature 
in the cortical structure which is responsible for formation of subjective surfaces. In the model they 
propose they study the evolution of a graph, and this is one of the main reasons why in this paper 
we focus on diffusion driven motion of graphs. 

The diffusion driven motion by curvature problem can be formulated also in the general setting 
of evolving arbitrary surfaces in Carnot groups (where it could have a wealth of applications in the 
study of minimal surfaces and isoperimetric problems). In this wider context, however, the notion 
of horizontal mean curvature is not yet known, because the presence of characteristic points cannot 
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be excluded. Although in this paper we do not study this more general setting we do provide 
a relation between mean curvature and heat kernels which could be used for that purpose (see 
Proposition [3]). 

In Section 2 we will start the proof of our result, by stating a preliminary technical result, which 
allows to express sub-Ricmannian quantities as limit of Riemannian ones. In particular we express 
the fundamental solution of the heat operator as limit of the corresponding Riemannian ones. This 
instrument is a generalization of the celebrated lifting method of Rothschild and Stein [47] , and in 
the time independent case has been developed in [14] . The adaptation to the heat operator provides 
a powerful and general instrument for transfer results from the Riemannian to the sub- Riemannian 
setting. The proof of the technical result can be found in the forthcoming paper |10j . 

In Section 3 we introduce the diffusion driven motion by mean curvature, and we restate the 
definition of the mean curvature operator in the setting of non linear semi-groups. In particular we 
recognize that the curvature operator A defined in fjl .5[) is non dissipative (see definition [2] below). 
According to a Theorem of Crandall and Liggett \17\ Theorem I, page 266] this allows to give a 
weak definition of solution of (jl.5p in terms of a non linear semi- group: 

Theorem. If A is non dissipative operator on a Banach space B, then for all f € B the limit 

(1.6) M{t)(f):= lim (I + \A)~> f 

j— yoo,Xj— H 

exists locally uniformly in t. 

The main result of section 3 is the proof that this notion of weak semi-group solution coincides 
with the viscosity one, given in Definition [TJ 

Proposition 1. For every continuous and periodic function f , the weak semi-group solution 
M{t)(f) is the viscosity solution u of motion by curvature, with initial datum f. 

Note that a non linear semi-group provides us with a discrete approximation of this motion, but 
the approximation is non-linear, while we were looking for a linear discrete approximation. 

In Section 4 we give a definition of the evolution operator H analogous to the one defined in 
[27j . However we are able to provide a substantial simplification, since we work only with graphs 
(and not with more general level sets), and we evolve only the subgraph of the given initial datum. 
Calling xs the characteristic function of its subgraph, we evolve it with the sub-Riemannian heat 
flow inGxR: 

(1.7) d t w=J2 X f w w(x,0) = X s(x). 

d(i)=l 

The new function H(t)f is defined requiring that its sub-graph is {x : w(x,t) < 1/2}. 

Note that the right-hand side of the heat equation (jl.7p is a sum of squares of m + 1 vector 
fields in R n+1 , and consequently it is only degenerate parabolic. However, in view of the work of 
Hormander [32], Rothschild and Stein [17] (see also Jerison and Sanchez-Calle [36]) such operator 
is hypoelliptic and admits a fundamental solution T with good estimates (see Proposition [5] and 
(|2.6p ). We conclude Section 4 with Theorem [2l which contains many properties of the operator H 
and in particular 

Proposition 2. For every continuous periodic function f and g 

(1.8) \\H(t)f - H(t)g\\ < \\f-g\\ forallt>0. 



HEAT KERNELS AND MEAN CURVATURE FLOW OF GRAPHS 



5 



Section 5 contains the main step and the geometric core of the proof: we will prove that the 
operator H locally approximates motion by horizontal mean curvature. The evolution of the char- 
acteristic function of a given set S is expressed through the heat kernel T as 



Following Evans [27], we need to perform a point- wise asymptotic expansion of w for t near 0, 
assuming that the boundary of S is sufficiently smooth. In the Euclidean case the expression of 
r is explicit, and the proof in \27\ Theorem 4.1] is a direct computation. In contrast, in our sub- 
Riemannian setting the heat kernel does not have in general explicit expression and in order to 
extend Evans' result we are obliged to use a deeper method, based on geometrical properties of the 
space. Our approach uses in full strength the properties of the sub-Riemannian heat kernel and its 
Riemannian approximation. In this extension, the results stated in Section 2 play a crucial role. 
Proceeding in this fashion we then establish some new formula, emphasizing a link between the 
heat kernel and the mean curvature, which are of independent interest. 

These new formula are closely related to the notion of heat content developed by De Giorgi [23], 
|24| in the Euclidean case, and recently extended to the sub-Riemannian case in [7]. In these works 
the integral of the function w is considered under low regularity on the set S, and it is proved that 
its derivatives tend to the perimeter. Higher terms expansions of the integral of w in the Euclidean 
case have been established in [3] where it has been proved that next term in the expansion depends 
on the mean curvature. 

Here we establish instead a pointwise estimate of the function w involving the curvature of OS 
Proposition 3. If S is the subgraph of a smooth function f and x £ dS, one has 



Note that this expression can be considered as a weak definition of horizontal mean curvature, 
which does not rely on differentiability of the surface or on other geometric considerations. The same 
argument provided in the proof of Proposition[3]yields an analogue result for arbitrary hypersurfaces 
of a Riemannian manifold (see Corollary [1]) . 

Remark 1. See the proof of p~[ Theorem 3.3] and |27] Theorem 4.1] for similar results in the 
Euclidean setting. In both papers the explicit form of the Euclidean heat kernel plays a key role. 

From the latter we infer: 

Proposition 4. For every continuous periodic function f one has 



uniformly as t — > 0. 

This approximation result allows to apply a general theorem of Brezis and Pazy [8] , which ensures 
the following: 

Theorem. If A is an m— dissipative operator {H(t)} t >o satisfying (|1.8p and 




2 Jn 

where II is the intrinsic tangent plane to S at the point x, as defined in the statement of Lemma\4\ 




t-\l-H{t)){f){x)^-Af{x), 



(1.9) 
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for every f and A > 0, then for every f £ D(A) and t > then one has 

(t\ j 

(1.10) M(t)f = lim H\ — ] /, uniformly for t in compact sets 



On the other side, since we have already proved the relation between the weak semi-group and 
the viscosity solutions of motion by curvature, we can finally deduce our main approximation result 
for the curvature flow: 

Theorem 1. If f is a continuous and periodic function, and u is the unique viscosity solution to 
(jl.5p . with initial datum f , then 

( t\ j 

u(x,t) = lim HI — I /, uniformly for x £ G and t in compact sets . 
j->oo \jj 

2. The sub-Riemannian structure as limit of its Riemannian approximation 

In this section we recall that the sub-Riemannian structure (G,go) in a Carnot group can be 
interpreted as a degenerate limit (in the Gromov-Hausdorff sense) of Riemannian spaces (G,g e ) as 
e — > 0, [30]. In [14] . Citti and Manfredini studied the relation between the sub-Riemannian Laplace 
operator and its Riemannian approximation. In particular they established bounds, uniform in the 
approximating parameter e, for the fundamental solution of the g e — Laplace operators. Here we 
recall a recent extension of these bounds to the degenerate parabolic case by Manfredini and the 
first two named authors |10j . 

We will work in the Carnot groups G and G = G x R, whose Lie algebras we denote by Q 
and by Q. On G one can define a semi-norm \x\^ = \(x, x n+ i)\^ := \x\& + a distance 

d(x,y) = and dilations 5\(x) = (6\(x), Xx n+ i) analogous to the same objects defined G. 

For each e > set 

X l = X\, ■ ■ ■ , = X m , = eX m+ i, ■ ■ ■ X^ = eX n , and X n+1 = X n+ i 

and V e u = (Xf )i=i,- , n - We extend both go and go to Riemannian metrics g e and g e on Q and Q by 
requesting that the vectors X\ are an orthonormal family. The degree of these new vector fields is 
d e (Xf) = 1 for e > and i = 1, ...,n + 1. The Carnot-Caratheodory distances (see [44]) associated 
to these vector fields can be proved to be equivalent to the pseudo-distances d e (x,y) := 
and d t (x,y) := associated to the non-homogeneous norms. For every v £ Q and v n+ \ G M, 

m n I, 

(2.1) |u|e=5^N+ Yl min (— >N lMi) ) and \{v,Vn+l)\e = Me + K+l|- 

i=l i=m+l 

It is then clear that for every e > one has d\(x, y) < d e (x, y) < d(x, y). Analogous considerations 
hold for d and d e . It can be shown (see |30] and references therein) that (G,d e ) — > (G,d) in the 
Gromov-Hausdorff topology. 

2.1. Riemannian mean curvature of a graph. As it is well known, the mean curvature of the 
graph of a regular function u t in the metric g e can be expressed in the form: 

1 n f 

h e = -= = = V a\AV e u)XfX* jUe with at JO = 8 k - 
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From this expression it immediately follows that h e — > ho as e — > 0. 

Also note that classical results establish existence, uniqueness and regularity for solutions to the 
elliptic equation 

n 

(2.2) u e - A a|,(V^)X|X> e = / 

for any regular domain of G, which can be identified with M. n (see |29j). 

2.2. Riemannian and sub-Riemannian heat kernels estimates. For e > 0, we denote by 
A e = Yld(i)=i X t the Laplace operator in G defined in terms of the metric g t . The subelliptic heat 
operator is defined in (jl.7p as a sum of squares of vector fields 

L = (d t - A )u 

while its Riemannian counterpart is 

L e = (d t -A e )u=(d t - X ?~£ 2 E X ^ u 

d(i)=l d(i)>l 

Both operators have fundamental solutions, which we call T (or occasionally To) and T e respec- 
tively. Note that the fundamental solution is defined on (G x K x M + ) 2 . However due to the Lie 
group invariance, we can fix the fundamental solution with pole in (0, 0), and the value in any other 
point will be obtained by translation invariance, i.e. for e > 

(2.3) F e ((y, t), (x, r)) = ^((x^y, t - r), (0, 0)). 

For simplicity we will drop the (0, 0) term in the argument of T and consider the fundamental 
solution as a function defined on M n+1 x IR + . Let us recall some well known properties of T, (see 
references in [36] ) 

Proposition 5. (1) For all x E GxK w have 

^t) = w ^T(5 t (x),l) 

(2) There exist constants C%, C2 such that, for all x G GxR and t > 0, we have 

(2.4) HQ+D/ 2 e -<* ^ < r(x, t) < t- W+ 1 )/ 2 e- c > ^ 

and 

iX'dfTfat)] < C 1 H«+ 1 )/ 2 -^- a e- c 2 ~ 1J2 ( £ > )/', 

where we have let I = i/), X 1 = X^.-X^ , d(I) = Y?h=i an( ^ Q = Yli d(i)dim(V l ) 

is the homogeneous dimension ofG. 

(3) For every t > 0, 

(2.5) / T(x,t)dx = l. 

JGxR 

This result implies in particular that for every fixed e > Gaussian estimates of the heat kernel 
T e hold with constants C\, C2 depending on e. In this paper we will need to use of uniform estimates 
of the heat kernels T e in terms of the limit kernel T (see |10j): 
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Proposition 6. There exist C\,Ci > such that for every e > one has 
(2.6) C7f 1 | J Be(0,\/t)r 1 exp(-C 2 4(x,0) 2 A) < F e (x,t) < Ci|fl e (0, V~t)\~ l ex^-C^ 1 d e {x , 0) 2 /t) 

and \{X e yd?T e (x,t)\ < C\\B e {0, v^T 1 *"^ 1-0 ex^-C^d^x, 0) 2 /t), 

where we have let I = (X 6 ) 1 = X^.-.X^ and d e (I) = X^L=i ^ e ■ Moreover, as e — > 

one /ias 

x/apr 6 -> x 7 af r 

uniformly on compact sets. 

The proof of this result can be found in [10] and is directly inspired to the arguments introduced 
in [TJ], where the time independent case is studied, and a general procedure (reminiscent of the 
lifting technique in [17]) is introduced for lifting the operator to a higher dimensional space. 



3. Nonlinear semigroups and horizontal mean curvature flow of graphs 

In this section we will restate the definition of mean curvature operator and viscosity solution of 
mean curvature flow in the setting of non linear semi-group theory. 

Definition 2. Let B be a Banach space, and A a B— valued nonlinear operator with domain 
D(A) C B. We say that —A is m— dissipative if 

• R(I + \A)=B for every A > 0, 

• its resolvent J\ = (I + A^4) _1 is a single- valued contraction 



The curvature flow will be defined for regular functions, periodic in the following sense: 

Definition 3. Let G be a Carnot group with a set of generators {e±, ...,e n } closed under Lie 
brackets, and denote by Q C G a fundamental domain, i.e. an open bounded set such that 
CiQnQ = for all % = 1, n and G = (J e" 1 ...e" n Q, where the union is taken among all n— tuples 
of integers a = (a%, ...,a n ). A function u : G — > ]R is periodic with respect to Q and e±, ...,e n if 
u(eix) = u(x) for all i = 1, ...,n. 

From now on B will denote the space of continuous periodic functions in a Carnot group G, and 
(B, || • ||) will be the Banach space obtained by endowing B with the sup norm || • ||. 

Definition 4. We say that u S B belongs to the domain of A if there exists / £ B and A > such 
that u is a weak solution (in the sense of Definition [5] below) of 

m , XuXu \ 

(3.1) u- aijiVoujXiXju = f, where aij(V u) := ISij - ^ . 2 J 

in G. In this case we will write 

(/ + \A)u = f. 

Definition 5. A continuous function u : G — > R is a weak sub-solution (resp. a super-solution) of 



(3.2) v-X^-J^y^f, 
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in G if for every x £ G and smooth <j) : G — > R such that u — 4> has a maximum (resp. a minimum) 
at x one must have 

u(x) - A ± - X ^ly { ^ y*XM*) < (resp. >) /(*). 

Solutions are functions which are simultaneously super-solutions and a sub-solutions. 

Proposition 7. The operator — A introduced in Definition^ is m—dissipative 

Proof. We first prove that \\u — v\\ < \\J\u — J\v\\. For A > let / = J\u = u + \Au and 
w = J\v = v + XAv, in the weak sense of Definition HI For ji > set u M and the sup-convolutions 
defined in 1521 : 



/ 1 \2r! / 1 \2r! 

u^ix) = sup ( u(y) d{x, y) ) , t> M := inf I v{y) + -d(x, y) ) 



In view of Wang's results in [52] one has that is a super-solution of (13. II) and is semiconvex. For 
each a,/i > we can then apply Jensen's lemma (see [37] or [161 Theorem 3.2]) to the semi-convex 
function in G x G, 

^W(z,y) := « M 0) - Vfj,(y) - ^\x - y\ 2 s 

and obtain a sequence of points Pj' a = (Xj' a ,yj' a ) —> p^' a = (x^' a , y^' a ) with p^' a a maximum 
point for the function which in turn is twice differentiable at pj' a and satisfies 

\V E <l>»,a(P? a )\ = (!) and - 7 J < D 2 E ^ a {pf a ) < o(l)I as j -> oo, 

where V# and Z?|, denote respectively the Euclidean gradient and Hessian, for some choice of 7 > 0. 
This implies that 

lim |*u"(xf) - XvM> a )\ = and lim XV(xf a ) - X 2 Vfl (y^ a ) < 0. 
j— too J J 

Since by definition u and v are solution of the PDE's f = u + XAu and w = v + \Av respectively, 
so that and satisfy 

m m 

- f < A ^(Vo'U^XjX^ at a^'" an d ^ - w > A ^ a ij (Vo^)^Q^j'iV at ?/^' Q . 

i,j=l ij'=l 

Subtracting the first inequality from the second and using the following fact proved in [161 Lemma 
3.2] 

lim a\x^' a — y^ ,a \ = and lim sup 4>a a (x, y) = sup(u(x) — v(x)), 

one concludes that ||u — u|| < ||/ — u>||, which is the first part in the definition of m—dissipative 
operator. 

We now establish that R(I + XA) = B. Let / € B D C°°(G) and for each e > consider weak 
solutions u t of the approximating elliptic PDE ()2.2[) . A simple variant of the argument above shows 
that for every e > one has \\u e \\ < ||/|| and, in view of the fact that the vector fields X[s are 
left-invariant, that 

Vy £ G sup \u t (yx) - u(y)\ < sup \f(yx) - f(y)\. 
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The latter proves that {u e } e is equi-continuous and equi-bounded. For a subsequence one has 
u ek — > u G B as efc — > oo and the argument in [9) Theorem 5.2] shows that u is a weak solution of 
(|3.ip . In view of the comparison principle we deduce that the range of / + XA is closed, hence it is 
the whole of B. 

□ 



The proof of the following proposition is exactly as the proof of the first part of [271 Theorem 
2.4]. 

Proposition 8. The domain of A is dense in B. 

Propositions [7] and [8] allow us to invoke the previously recalled Generation Theorem of Crandall 
and Liggett (see 111.6)) above), and ensure that —A generates a non linear semi-group M(t)(f). 

Proof of Proposition [1] 

Let (f> : G x (0, oo) — > R be a smooth function such that u — (f> has a maximum point at 
(#o>*o) G G x (0, oo). We can always assume that the maximum is strict adding a power of the 
gauge distance as in [9]. We set 

/ 1 \ r k k + I s 

u s (x,t) = ( (I + -A)~ fc /)(x) for all integers s, k and i G 



s 



s s 



and note that for a fixed t > 0, as s — > oo one has fe — )• oo and fc/s — > i. By Crandall and Liggett's 
convergence result one has that u s — > u uniformly on compact sets. As a consequence one can find 
a sequence (x s ,t s ) G G x (0, oo) such that: 

(i) (x s ,t s ) -> (x ,io); 

(ii) For fixed i s > the function (u s — </))(■, t s ) has a maximum at x s ; 

(iii) {u s - (/)(x s ,t s ) + jr>{u s - f){x,t) for all (x,t) G G x (0,oo). 

Following [27) Theorem 2.5] for each s we let k s be such that t s G [— , ^+i) and in particular 
t s -i£ [^fi, %). This choice yields 



u s (x,t s ) = (I+-Ay 1 u s (;t s --)(x) 
s s 



in the weak sense of Definition By virtue of (ii) and (iii) above one has 

c/)(x s ,t s ) - 4>{x s ,t s - \) ^ / Xi(j>Xj<i> \ . 1 



^/ \ 1 



i 

S jj = l 

As s — > oo, the latter completes the proof that u is a sub-solution. In a similar fashion one can 
prove that u is also a super-solution and conclude the proof. 

4. A DISCRETE TIME-STEP OPERATOR BASED ON THE HEAT FLOW 

4.1. Flow of subgraphs of functions belonging to B. Let (B,\\ ■ ||) be the Banach space 
introduced in Definition UJ If / G B, we call subgraph of B the subset of G x R defined as 

(4.1) S = {x = (x, x n+1 ) : x n+l < f(x)} C G x R. 
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We give here the definition of heat flow of a subgraph S defined in (14. ID . as an extension of the 
analogous notion given by Evans in the Euclidean case |27] , For all t > denote w(x, t) the solution 
of the Cauchy problem 

Aw = dtw 
«>(•, 0) = xs 
and define 

(4.2) U(t)(S) = {x G G x R : w(x,t) > ^} 

In order to prove that this functional induces a flow on B, we need to show that if the subgraph 
of a periodic, continuous function / : G — > R is evolved through the flow H(t) then the new set is 
also the subgraph of a periodic continuous function. 

Lemma 1. Let / G B. For any t > there exists a function g G B such that 

H(i){x n+1 < f(x)} = {x n+ i < g{x)}. 
Proof. Let (y,y n +i) & d{H{t){x n+ i < f(x)} so that 

-= T(x~ 1 y,t)dx = / / r(x~ 1 y,y n+ i - x n+ i,t)dx n+1 dx. 

1 Jix„+i<f(x)\ JgJ-oo 



2 

If we denote 

r/(*) 



(4.3) F f (y,y n+1 ) := / / T(x y n+ i - z n+ i, t)dxdx n+ i. 

J <Q J — oo 

differentiating along the variable y n .+i yields 

d y n +i F f(y>yn+i) = / r(x _1 2/,2/ n+ i - f(x),t)dx < 0. 



The implicit function theorem then implies that in a small neighborhood of every point z G G, 
there exists a continuous function g : B(z,r) C G — > M such that Ff(y,g(y)) = 1/2 for all 
y G -B(z, r). On the other hand, by the strict monotonicity of the function Ff as a function of its 
last variable, we immediately deduce that the function g is defined on the whole group G with real 
values. 

□ 

This lemma allows to give the following definition of flow of functions in B: 

Definition 6. If / G B we define the operator 

H(t):B^B, H(t)(f)=g, 

where g is the unique function satisfying 

U(t){x n+ i < f(x)} = {x n+ i < g(x)}. 

Thanks to the explicit expression of H we can establish comparison principle properties of the 
flow H. 

Theorem 2. For each t > the flow Hit) :B^B just defined has the following properties 

(1) If f <g thenH(t)f <H{t)g 

(2) If C is a real constant, then H(t)(f + C) = H(t)f + C 
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(3) H{t) is a contraction on B 

\\H(t)f-H(t)g\\<\\f-g\\ 

Proof. Let us show assertion (1). If /, g £ B, and / > g we can use Lemma [H and represent the 
function F g defined in (|4.3p as 

r rg( x ) 

F g (y,y n +i) = F f (y,y n+1 ) +1 l T(x~ y, y n+ i - x n+1 , t)dxdx n+1 . 

J& J f{x) 

Since from the definition of Ff and H(t)f we have 

F f (x,H(t)f(x)) = F g (x,H(t)g(x)) = 1/2, 

we immediately get 

F f (x,H(t)f(x)) = F g (x,H(t)g(x)) > F f (x,H(t)g(x)). 
By the strict monotonicity of Ff in its last variable we deduce: 

H(t)f < H(t)g. 
In order to prove assertion (2) we note that 

r rf{x) 

Ff(y,y n+ x)= / T(x y,y n+1 - x n+1 ,t)dxdx n+ i = 

J& J —CO 

with the change of variable z n+ \ = x n+ \ + C, 

"f{x)+C 

T(x l y, y n+ i + C - z n+1 ,t)dxdx n+1 = F f+C (y, y n+ i + C) 
For y n+ i = H(t)f(y) we get 

F f+ c(y,H(t)f(y) + C) = 1/2, 

which implies by the uniqueness that H(t){f)(y) + C = H(t)(f + C)(y). 
Finally we prove assertion (3). Let us choose x such that 

H(t)f(x) - H(t)g(x) > \\H(t)f - H(t)g\\ - e 

and call 

H(t)f(x) = n H(t)g(x) = X. 

By assertion (2) we have 

H(t)f(x) - H(t)(g - A + ii - e)(x) > 
By assertion (1) this implies that there exists a point y such that 

f(y)-(g(y)-\ + pt-e)>Q, 

which by definition of A and [i implies 

\\H(t)f - H(t)g\\ < H(t)f(x) - H(t)g(x) + e < f(y) - g(y) + 2e < | \f - g\ \ + 2e 
Assertion (3) follows. □ 
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5. SUB-RlEMANNIAN DIFFUSION OF SMOOTH NON-CHARACTERISTIC GRAPHS 

In this section we prove the geometrical core of the paper, which is the equality of the normal 
speed of a graph evolving by mean curvature and the normal speed of the same graph evolving under 
the heat flow defined in (14, 2\i . According to the notation introduced in section 2, for / € B we will 
denote by do~\ the surface measure element on dS = {(x, f(x)) 6 GxM|x G G} induced by the metric 

g e , with e = 1 and by dao the horizontal perimeter measure dao = V 1 ~ t ~j v °" ]I ( fo- 1- T/he outer hori- 
zontal unit normal v§ to S is expressed in the frame X\, ...,X n+ i as v$ = (— Vo/, l)/yl + |Vo/| 2 . 

Lemma 2. If / : G — > K is a smooth function, with sub-graph S, then for any smooth functions 
u,v £ B and for every i = 1, m one has 

v(x)Xiu(x)dx = — / u(x)Xju(x)dx + / u(x)u(cc)^do"o(a;) 
Proof. Integration by parts yields 

v(x)Xiu(x)dx = — ( u(x)Xiv(x)dx + / u(x)i>(x) — = =do~i(x) 
Js Jds -y/1 + IVimI 2 

u(x)Xiv(x)dx + / ii(x)-u(x) — 1 =rf<7o(x). 
5 ./as Vl + |V /| 2 

□ 

Lemma 3. Let B be the space of periodic functions defined in Definition [H If f £ B let S its 
sub-graph (defined in (|4.ip . and let x = (x,x n+ i) = (x, /(x)) be a fixed point in <9,S. Given L > 
for every u; £ Lip(G x M) with Lipschitz norm bounded by L, one has 

/ T(x- l y,T)w(y)da (y) = ^ [ T(z,l) da (z) + 0(1) 
J as v r Jn 

as r — 7- uniformly for x £ dS. Here II is the intrinsic tangent plane y n +i = n(y) := /(x) + 

Z)d(i)=i x if( x )(Vi ~ x i) to dS" at x. 

Proof. Since dS is a graph, the integration variables are of the form y = (y, f{y)). Note that 

(i- 1 y) = (x- 1 yj(y)-f(x)) = (x- 1 y, ^ Xif(x)(y - x\ + 0(\y- l x\ 2 )). 

d(i)=l 

Consequently, setting z = S T -i/ 2 (x~ 1 y) one obtains 

(5.1) z = (z,z n+1 ) = (5 T - 1/2 (x- 1 y),r- 1 / 2 (f(y) - /(x))) 

= (K^{x- l y), Xifix^-y^y^ + O^z] 2 ) = (z,U(z) +0{V^\z\ 2 )). 
d(i)=l 

In view of the latter and using the scaling properties of the heat kernel one has 

(5.2) T(x^y, r) = r^ 2 " 1 / 2 !^, 1) = T ~ Q ^ 2 T((z, U(z) + 0(^\z\ 2 ), 1) 



r- Q/2 ~ 1/2 (r((z,n(z)),i) +J R(z)) 
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where R is a remainder satisfying f G Rdz = 0(^\z\ 2 ). Consider the surface integral 

(5.3) / T(x~ 1 y,T)w(y)da (y) 
JdS 

= w(x) T(x~ l y,T)da (y) + 0[ \x~ 1 y\T(x~ 1 y, r) da Q (y) ) 

JdS WgxM 7 

= w(£) / Tix-iyJiy) - x n+1 ,r)y/l + \V f(y)\ 2 dy + 
JG 

+ 0(J \(x- 1 y, f(y) - x n+l )\ T{x~ l y, f(y) - x n+1 , r) y/l + |V /(y)| 2 dy) 
Applying the change of variable z = 5 ^(x~ l y) we have dy = T®' 2 dz and obtain 

(5.4) / r(x-V f(y) - x n+1 , T )y/1 + \V f(y) \ 2 dy 

Jg 

= ^J ( r ((^ n (^)) 5 1) + ^l + \V f(x + 5^(z))\ 2 dz 

= ^J g (r((z,U(z)), 1) + R(z)j v 7 ! + \VoWdz +0(1) 

= 4= / r(5, 1)^0(^ + 0(1). 

where in the last line we have used the approximation V/(x + 5 ^p(z)) = V/(x) + 0(\A"|^|) and 
we have denoted z = (z,Tl(z)). The proof follows immediately from the latter and from (|5.3H . □ 

We are now ready to establish the key geometric identity needed in the proof of Proposition [3] 
and eventually to prove Theorem [TJ 

Remark 2. The previous lemma indicates that it would be plausible to have an identity of the form 

/ \VT(x- 1 y,r)\w(y)da (y) = ^^- + O(^ r ) 
JdS T V T 

However in the next lemma we show that for integrands with a special structure one can improve 
on such an identity an obtain a better decay rate. 

Lemma 4. Consider / € B a C 3 function and suppose that the horizontal mean curvature ho of 
its graph is bounded by a positive constant L. Let S denote its sub- graph, defined in (|4.1|) . and 
consider x = (x,x n+ i) = (x, f(x)) a fixed point in dS. Then one has 

(5.5) / (V r(x- 1 y,r),u (y))- go dao(y) = -^ [ Y(x~ x z, l)da {z) + O(l) 
Jds 2 V T Jn 

as r — > 0, and the convergence is uniform with respect to the variable x £ dS, but depends on the 
parameter L. Here n is the intrinsic tangent plane y n +\ = n(z) := f(x)+Yld(i)=\ ^if{ x )iVi~ x i) *° 
S at the point x, dao denotes the horizontal perimeter measure of n associated to X\, X m ,d n +\. 
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Remark 3. Note that if x £ dS then there exists a = a(S, x) > such that 

(5.6) / (V T(x- 1 y,T),u (y)) §0 da (y) = O(e- a / T ). 

JdS 

In particular the left-hand side is uniformly bounded by 0(r -1 / 2 ) for any point ieCxI, 

Proof. The proof proceeds in two steps. We will start considering the extremely simplified case 
in which the group G is the Euclidean space and the horizonal tangent plane coincides with the 
whole tangent plane, so that the sub-Riemannian metric reduces to a Riemannian one, and m = n. 
For reader convenience we will keep also in this case the same notations used in the rest of the 
paper. In the first step the metric will be constant, in the second step we consider an arbitrary 
Riemannian metric, while the subriemannian case will be studied in the third step as limit of the 
Riemannian ones. 

Step 1 Let us first assume that G = 1" with the Euclidean group structure and that we have 
a constant coefficient Riemannian metric go defined on the whole tangent space of G x E. Up to 
a change of variables we can assume that x is the origin, that (X\, ■ ■ ■ X n ) = (d yi , ■ ■ ■ ,d Vn ) is an 
orthonormal frame for Tq8S and X n+ \ coincides with the outer unit normal vq to S at 0. In these 
coordinates the matrix associated to the metric go is the identity, the heat operator (defined in 
Section 2) simply reduces to the standard Euclidean one dt — Xf. 

There exists a neighborhood U C Gxlof the origin where U(~)dS can be represented as a graph 
of a smooth function over U n TqOS = U n {y n +\ = 0}. We will still denote by / this function and 
by (y, f(y) its graph. The normal is expressed fo(y)dao(y) = (— Vo/(y), and, by the choice of 
coordinates, Vo/(0) = 0. It follows that 

{V Q T(5;- l y:,T),v Q {y))g Q doo(y) = [ (Vol^y, r), v Q {y))~ go da {y) + 0( e ~ a /*) 
dS Jasnu 

(where a > is a suitable constant) 

<(y,/(y)),(-Vo/(i/),D)/ r(yi f{y)t)dy + 0{e - a/t) 

unT ds zt 

where (•,•)/ denotes the Euclidean scalar product. Extending in a periodic fashion / to all of T$dS 
and identifying the latter with W 1 one obtains 

m-V 'J^ T(yJ(y),t)dy + 0(e-^) 

The integrand function can be expanded as follows: 

f(y) - yV f(y) = ^WO) + o(\ y \ 3 ) 

Applying the change of variable y = Viz, and the scaling property of the heat kernel in R n+1 , we 
arrive at 

/ (V T(y,r),u (y)) §0 dao(y) = ^/P-( f ^ T(z, 0, 1) dz + 0(V~t)) 
Jds Wt \Js. n 1 J 

= &^L I r(z,0,l)da (z) + O(l) = -^ [ T(z,l)da (z) + O(l), 
Wt Jn Z V T Jn 
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here we have used the fact that in this system of coordinates ho(0) = —daf(0). This concludes the 
proof in the constant metric case. We explicitly remark that in the setting of step 1 the hyperplane 
n coincides with the tangent plane to S at x. 

Step 2 Let us now assume that the metric go = E T E is an arbitrary Riemannian metric on 
invariant with respect to a Lie group G. We can define an orthonormal frame in R n+1 by 
setting Xi = (E~ 1 )fjdj, with i,j = l,...,n + 1. The associated heat operator now contains first 
order terms 

71+1 

L = d t - X? = dt - ~g kj d% - E^djiE-^ldk = d t - ~g^d% - ~a k d k 

i=l 

where {g^} denotes the inverse of g, and d k = EJ^dj^E -1 )^. As in step 1 we can assume without 
loss of generality that Vo/(5) = 0. Consequently — a k (x) = and the associated constant coefficient 
operator, obtained evaluating the coefficients at the point x, reduces to: 

n+l 

L S = dt-Y^Xl^ d t - g kj {x)d%, 
i=i 

where we have denoted Xg^ the vector fields with coefficients evaluated at the point x. We will 
denote by = (X^i, • • • X^^+i) the gradient along such frozen coefficients vector fields. We will 
also denote by V and the heat kernels of L and Lj. Finally we denote by ho and ho,x the mean 
curvatures of the graph with respect to the metric go and with respect to the frozen metric g(x). 
Note that at the point x one has hx(x) = h(x). 

We want to emphasize that in this step of the proof we are never going to use the invariance 
properties of the heat operator L and its fundamental solution T. Because of this, the group 
structure itself becomes irrelevant and the fundamental solutions will be considered as a function 
defined on (M n+1 x M+) 2 , thus requiring the notation T((y,t), (x,t)) (see also discussion at the 
beginning of Section 2 and (j2.3[) ) . 

Using the Levi's parametrix method (see for instance [32], |28j). it is known that one has 

(5-7) T((y, t), % r)) = T,((y, t), (x, r)) + o( _ 1 e' C ^) . 

More precisely 

(5.8) r((y, t), (x, r)) = r 5 ((y, t), (x, r)) + J((y, t), (x, r)) 

where 

J((y,t),(x,r))= [ r s ((y,t),(z, S ))Z 1 ((z : s),(x,T))dzd S + o(- i_^ e -^^#) 

and 

Z 1 {{y,t),(x,r))=Y,{9 lJ {y)-~9 l3 (x))d l] T, i {(y 

ij i 

so that, for some constant C = C(G,g) > 0, 

,r((y,i),(x,r)) 



\Z l ((y,t),(x,r))\<C- 
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Clearly one also has 

(5.9) V r((y, t), (x, r)) = V r £ ((y, t), (x, r)) + V J((y, t), (x, r)). 
Let us consider the integral 

(5.10) / (V J((y,t),(x,r)),i/ (y))g ^o(y) = 



( / VT B ((y,t),(z,s))Z 1 ((z,s),(x,T))dzds \, i/ Q (y))g o da (y) 

dS WR"+ 1 x[r,t] / 



(V r i ((y,t), (2,s)),i/ (y))g dcro(y) )Zi((z,s), (x,r))dzds. 

l n + 1 x[r,t] y JdS 

In view of Remark [3] there exists C = C(G, g, S) > such that 



(V J((y, t), {x, r)), My))g dao(y) 

as 



< C 



\M(z,sUx,r))\ d2ds + 0{e _ a/ ^ T)) 

"+ 1 x[r,t] Vt— S 

< C f * = / r,((y,t),(5, S ))^ + 0(e- a /^)) = 

•/ t v(t - s)(s - T)) JR"+i 

Wr V(i-«)(s-r)y Wo Jr(l-r)J' 



o 



V / (t-s)( S -T)J \J y/HT 

where, in the last line, we have used the change of variable r = (s — r)/(t — r). 
Consequently, 

(V r((y,t), (x,T)),v (y))g da (y) = 

dS 

(V r £ ((y, t), (x, r)), is (y)}g da(y)+ 

dS 

+ jj V * J ^i ')> & T ))' ^(y))so^o(y) + ° ^ s ( t _t)n/2 e " C2 ^ d<7 ^)) • 

Denoting by da g (%\ the surface measure corresponding to the frozen metric gx) and using the 
statement in Step 1, one has 

(5.11) / {V T((y,t),(x,r)),u (y)y go da (y) = --^p= I T B ((z, 1), (5, 0))dtr g{i) (y) + O(l) 

(using Euclidean rescaling) = — / r £ ((jr,t - r), (x, 0))dcr fl m(y) + O(l) = 

z Jn 

(by estimate (|5.7p ) 

(5.12) - / T((y, i - r), (x, 0))d<r(y) + O(l) 

We explicitly remark that also in the setting of step 2 the hyperplane n coincides with the tangent 
plane to S at x. 
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Step 3 Let us now assume that T is a sub-Riemannian heat kernel corresponding to a sub- 
Riemannian metric go in the Carnot group G. From Section [2] one has a sequence of left-invariant 
Riemannian metrics g e in G such that (G, d t ) — > (G, do) in the Gromov-Hausdorff topology and 
corresponding sequence of Riemannian heat kernels T e satisfying the uniform estimates in Propo- 
sition El Applying step 2 at every level e > and using the fact that g t are left-invariant, one has 
the identities 

(5.13) f (V e T e (x- 1 y,T),u e (y)}gJ^(y) = -^4r- [ ^(^,^(^ + 0(1), asr ^0, 
JdS z JTidS 

where the bounds in 0(1) are uniform in e. If we denote by e the heat kernel corresponding to 
the frozen Riemannian metric g t (x) (no longer left-invariant with respect to G) then (|5,lip and 
([5TT2]) yield 
(5.14) 

/ (Vj: t {x- l y,T)My))~ gt doM = - f} ^ [ T ite ((z,l),(x,0))da ge{£) (z) + O(i), as r -> 0, 
JdS l V T JTzdS 

where the bounds in O(l) are uniform in e. 

In view of Proposition [6] and the dominated convergence theorem it follows that the left-hand 
side integrals 

/ {V e Y e {x- x y,T),u e {y))~ ge da e {y) 

JdS 

converge to the left-hand side of (|5.5p as e — > 0. 

Next we turn our attention to the right-hand side of (|5.5p . First we recall that at every point 
h e — > ho as e — )■ and da e — > dao as measures. Using the uniform bounds in Proposition [6] along 
with dominated convergence we deduce that 

lim / T e (y,t)da t (y) = T(y,t)da {y) 
t^JTodS JT dS 

(applying Lemma [3]) = — j= / F(z, l)dao(z) + O(l), 

V r Ju 

□ 

which completes the proof. 

Remark 4. The result in Lemma [5] seems to be new even in the Riemannian setting (Step 2 in 
the previous proof). Since we find that such extension may be of independent interest we state it 
explicitly. 

Corollary 1 . Consider a Riemannian manifold (JV, g) and a smooth embedded hypersurface M C N 
endowed with the induced metric. Denote by T be the heat kernel on (N, g), and for x £ N denote 
by the heat kernel in corresponding to the frozen metric g(x) viewed as a function defined 
on (TxN x M + ) 2 . We also consider da g , the induced volume element on M, the unit vector field v 
normal to M, and by h the mean curvature of M. For For every x 6 M and t > s > one has 

(5.15) / (V g T((x,t),(y,s)),u(y) G da g (y) = - H j-L f T s ((x, t), (z, 8 ))dffz{z) + O(l) 
Jm 2 Vt - s Jt £ m 

as t — > s uniformly for x, £ M. Here we have denoted by do~x the volume element on T X M induced 
by the metric g{x) on T X N. 
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The proof follows closely the argument in Step 2 of the previous Lemma HI As we already 
noted, in that argument we never used the group law structure, and Step 2 can be applied for any 
Riemannian metric and independently of the presence of a group structure. The result is local and 
we used the splitting of the space in G xl only to express the boundary of S as a graph. Since 
this can be always done, under suitable regularity assumptions, the result in Step 2 holds in any 
Riemannian metric. 



We can now prove the asymptotic expansion stated in the Proposition [3] in the introduction 

Proof of Proposition [3] Let p e be a smooth mollification of xs- For any interval (£1,^2) C M + , 
from the definition of heat kernel, 

(5.16) 0= lim/ 2 / (d t -C)T(^ 1 y,T)p e (y)dy= ( [ T(^ l y,T)p^y)dy) t2 

R-^oeJt! J{\x- 1 y\<R} V JgxR ' *i 



- lim [ [ (V r(5 1 y,T),v (y)} p € (y)da (y)dT 

+ r [ (V r{x- 1 y,T),V p e {y)}odydT, 
Jti JGxR 

Letting t%, e — > 0, one obtains 

(5.17) 0= / r(x- 1 y,t 2 ) X s(y)dy- lim / T(x" 1 y, ti)xs{y)dy 

~ i im / I (^o^(x~ 1 y,T),u (y))oxs{y){y)da {y)dT 

R-^°oJ0 J\\x-^v\=R\ 



t-2 r 

/ {V T(x~ 1 y,T)^ (y)} da(y)dT, 
J as 



Let us note that since x £ dS then 



lim/ r(x 1 y,t 1 )xs(y)dy = - 



Hence 

rt 2 



[ T(x 1 y,t 2 )xs{y)dy = i - lim f [ (V r(x 1 y,T),v (y)) xs(y)d<Jo(y)dT 

JGxM 1 R-+°°Jo J{ \x- x y\=m 



+ I 2 [ (V T(x- l y,T),v (y)) da(y)dT, 
JO JdS 



Next we show that 



(5.18) lim/ 2 /' (VolX^y, r), u (y)) X s(y)da (y)dr = 0- 

R-^°°J0 J{ {x-i-y^R} 



20 LUCA CAPOGNA, GIOVANNA CITTI, AND COSIMO SENNI GUIDOTTI MAGNANI 

To see this we recall that the dao perimeter of dB(0, R) in G x R is given by fg B r R \ dao = C<qB9 . 
From this and from the heat kernel estimates one has 



f (V r(x- 1 y,T),z. (y))o xs(y)da (y)dr < C G R~ 2 [ 
J{\x-^v\=m Jo 



-t 2 / R 2\{Q+2)/2 r2 

e~~dr, 

<{\x-iy\=R}~ Jo V T 



which implies (|5.18p . 

Applying Lemma E] one concludes 



f t2 (h (x) 



(5.19) / Y{x- 1 yM)xs{y)dy = ^- I [ I T(z,l)da (z) + 0(1) )dr 
Jgxr 2 Jo V 2 V T Jn J 

(5.20) = \ - h (x) [ T(z, l)da (z) [ % -^dr + 0{t 2 ) 

2 Jn Jo 2 V T 

(5.21) = \ - h (x)Vh [ T{z, l)da (z) + 0(t 2 ). 

2 Jn 

Theorem 3. Let f € B be a smooth function, and denote S its subgraph, defined in (I4.1h . If x £ 

x = (x,f(x)) a G M m+1 and t>0, denote 

q(t) = (x,f(x))exp( ^2 

d(i)=l 

Then we have 

q{t) G d{U(t)S) =► a i u o = ~ho{x)t + 0{t 3/2 ) ast^O, 

d(i)=l 

where ho(x) is the horizontal mean curvature of dS in x and vq = (z/q, Vq +1 ). 
Proof. From the definition of the heat flow of sets one has 



Y{q{t)- L y,t)dy 

considering the Taylor expansion of the integrand with respect to v we obtain 
= I T(x- l y,t)dy- a i [ X t \xT(x- 1 y,t)dy + o(t) 



d(i)=l 



/ F(x 1 y,t)dy- V" on I Xi\ x T(x 1 y,t)dy + o(t) 
Js ,777, Js 



d(i)=l 



moreover applying Lemma [2] we get 

= f T{x- l y,t)dy- V en [ r(x~ l y, t)v?(y)da Q (y) + o(t). 
Js ,77T, JdS 



d(i)=l 
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and applying Lemma [3] 



/ T(x- l y,t)dy- V / T(z,l)da (z) + O(t) 

Js d(i)=i vt Jn 

y/tho(x) + r(5, l)da (2) + O(i)) 



d(t)=l 

2 

(5.22) = 1 - -L jf r(z, l)d7 (5) («H)(5) + (a, ^o) 30 + <?(*>/*)) 

where we have applied Lemma [3] in the first integral, and Lemma 0] in the second, concluding 
the proof. □ 



Next we derive two important Corollaries from the previous theorem, which will be the main 
ingredients in the proof of Theorem [TJ 

Corollary 2. Choosing a = tfivQ, with f3 G M in the previous Theorem, we deduce that, if 
q(t) = (x)exp(tpvo) G d{H(t)S) theifl 

/3 = -ho(x) + 0(y/i) as t -)• 0, 

where ho(x) is the horizontal mean curvature of 55 in x 

Proposition 9. Xei f £ B be a smooth function and for t > denote by H(t)f its flow defined in 
Definition® For every x £ G one has 



(H(t)f)(x) - f{x) = -th Q (x, f(x))y/l + | V /(x)| 2 + o(t), 
where the convergence o(t)/t — > is uniform as t — > 0. 
Proof. Choosing a = ((H(t)f)(x') — f{x'))X n+ i we have 

(H(t)f)(x) - f{x) =< a, uq > go \/l + |Vo/(x)| 2 = -th (x, fix)) y/l + | V /(x)| 2 + 0(*) 

□ 

We can now conclude the proof of the main result of the paper, Theorem [H As in [27] the key 
technical tool in the proof is the non-linear version of Chernoff 's formula established by Brezis and 
Pazy (see Theorem II. 3D in the introduction. 

Proof of Theorem [TJ We only have to show (ll.9p for A = 1. To this end we set for t > and 



I + t-\l-H{t))\ f, and A t u:= U ^ } " . 



In view of Proposition [2] the operator —A t is m— dissipative, thus implying that for all y G G and 

t > 0, 

sup \u\yx) - u\x)\ < sup \fiyx) - f(x)\ 



x6G xeG 



and consequently that {"U*}te(o,i] is a family bounded and equi-continuous. 

Let <f) G C°°(G) such that u — <j> has a positive maximum at xq G G. We can always assume that 
the maximum is strict, adding a suitable power of the gauge distance, as for example in [9]. Since 



*Note the difference in sign with |27l Formuia (60), Theorem 4.1], which arises because of our choice of unit normal 
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u k —> u uniformly on compact sets then one can find a sequence of points x^ — > xq as k — > oo such 
that u tk — 4> has a positive maximum at x k . In view of Proposition [2] one has 

(H(t k )u tk )(x k ) - (H(t k )<f> tk )(x k ) < u tkk {x k ) - 4>(x k ), or equivalently A^<f>{x h ) < A^u tk {x k ). 

Since u l + A l x^ = f then 

(5.23) u*-( Xh ) + ^)-(ffa)^)^) < /(xfc) . 

Invoking Proposition [9] with <fi in place of /, one obtains 

Letting fe — )• oo we establish that u is a weak sub solution of (|3.2p with A = 1. In a similar fashion 
one can prove that u is a weak super-solution, concluding the proof. 
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